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The NMR relaxation rate l/Ti is studied for one-dimensional multicomponent spin-orbital 
systems. By combining the bosonization techniques and the exact solution of the SU(n) model, 
we evaluate the power-law exponent and the enhancement factor for 1/Ti at low temperatures. 
We discuss how the band splitting affects the relaxation rate, and find that l/Ti may be enhanced 
around the critical value of the band splitting. The crossover behavior in l/Ti around the critical 
point is discussed in terms of the low-frequency dynamical spin susceptibility. The effect of the 
hole-doping is also addressed. 
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§1. Introduction 

Quantum spin systems with orbital degeneracy have 
attracted considerable attention recently, since it was 
recognized that the orbital degrees of freedom give rise 
to a variety of interesting phenomena in correlated elec- 
tron systems. These hot topics have stimulated the- 
oretical studies on the orbitally degenerate Hubbard 
model and the related multicomponent spin-orbital mod- 
els. As a first step to investigate such orbital effects, the 
one-dimensional (ID) version of the Sli(4) spin-orbital 
modftl rh%fj ^Pp^ studied numericallyoB' and analyti- 
cally.B'cHj'O'LlB Such analyses have been|-|ais|0--e|5sl(pjded 
to the SU(2)xSU(2) spin-orbital modelMBBByM 
However, these studies have been mainly concerned with 
the phase diagram and/or the static properties, which 
naturally motivates us to study dynamical quantities for 
the above models. 

The NMR measurement may be a powerful probe to 
study low-energy spin dynamics in ID strnngl/y— ppci 
lated systems, as demppstLated theoreticall5iil§ll3'll3^t§ 
and experimentally££l^E2l'E3 In this paper, we study the 
NMR relaxation rate for ID multicomponent spin-orbital 
systems. For this purpose, we first give a multicompo- 
nent generalization of the formula for the NMR relax- 
ation rate 1/Ti in terms of the bosonization method. 
We then evaluate the power-law exponent and the en- 
hancement factor for 1/Ti by exploiting the integrable 
SU(n) spin-orbital model. It is discussed how the band 
splitting affects the relaxation rate. We further study the 
dynamical spin susceptibility to discuss the crossover be- 
havior in 1/Ti around the critical point. The effect of the 
hole-doping on the relaxation rate is also discussed. 

This paper is organized as follows. In §2, we obtain 
a generalized expression for the NMR relaxation rate, 
and then discuss the the effect of the band splitting on 
1/Ti in §3. We present the results for the dynamical spin 
susceptibility in §4 to discuss how the crossover behav- 
ior emerges when the band splitting is changed. Brief 



summary is given in §5. 

§2. Relaxation Rate for Spin-Orbital Systems 

2.1 Bosonization approach 

We first derive a generalized expression for the NMR 
relaxation rate 1/Ti for ID multicomponent electron sys- 
tems. To illustrate the procedure clearly, let us consider 
the degenerate Hubbard model with the band splitting, 
which possesses iV-fold orbitals and spins a =1, |, 



N 



h.c.) 



U 



tEE e 



^EEo 



(1) 



where c , is the creation operator for an electron with 
the orbital labeled by n (= 1, 2, • • • , N). The band split- 
ting A, which is assumed to be caused either by the 
crystal-field effect or by the Zeeman effect, gives rise to 
the energetically separated bands with the dispersion, 
e(fc) — —2tcoskao ± A where ao is the lattice spacing. 
We deal with the case that the lower (upper) band is 
p-fold (q-fold) degenerate, for which the factor 2 due to 
the spin sector is included. The number of electrons is 
denoted as n\''^ ("■i"'') for the lower (upper) band. Note 
that the system possesses SU(2A^) spin-odpital symme- 
try in the absence of the band splitting .Bl^ For con- 
venience, we use the number n ~ 2N and the index 
m — 1, 2, • • • , n{— 2N), which specifies the spin as well 
as the orbital degrees of freedom on equal footing. 

Passing to the continuum limit, we now introduce the 
right- and left-going fermion fields J£3' 

c™ ^ [e'^-^^P+M + e-''=-^V-m(a:)] , (2) 

where the Fermi momentum, kp represents either kp^ for 
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the lower band (to — 1, ■ ■ ■ p) ot k^p^ for the upper band 
(to = p + 1, • • • , n). These fermion fields are bosonized 
in terms of and its dual field 0™ , 



^rrn 



:exp(-i[r$,„ - Q„i]), 



(3) 



where the boson fields satisfy the commutation relation, 
[^m{x),Q^'{x')] = i(7r/2)(5„j„'Sgn(a;' - a;). The anti- 
tion relation between fermions with different 
is preserved by Majorana-fermion operator, 



comm 
specie 



Let US now bosonize the Hubbard interaction. We con- 
sider a metallic case, for which the charge, spin and or- 
bital excitations are all massless. In this case, we can 
neglect both of the backward scattering and the Umk- 
lapp scattering around the low-energy fixed point. This 
simplification reduces the system to an effective multi- 
component Tomonaga-Luttinger (TL) model in external 
fields, 



H 



E 



^ ^V,n [dx^mix)f + V„i [dxQm{x)f^ 



dx [dx^„,{x)] [dxKi'ix)] 



(4) 



where = vi for to — 1 , • • • , p and — for 
TO = p + 1, • • • , n are the Fermi velocities of elementary 
excitations. 

Note that the Fermi velocities in the lower and up- 
per bands are different, making the diagonalization of 
the Hamiltonian (4) a little bit complicated. Neverthe- 
less, we can perform the diagonalization via the canonical 
transformatioo. following the way done for a two-channel 
TL modelEil'E3 To this end, we introduce the new boson 
fields, <i>c, <&A and <i>cr,„ for the charge, orbital and spin 
sectors, respectively. For the lower band, the transfor- 
mation is given by 



1 , y/m . 

— cos a + y sm a 
py'n 



1 sin a 



pq 



In y 



P\fn 



cos a 



p+1 

E 



1 



„f^3 v/K-i)K-2) 

X \B(rn! - TO - i) - (to' - 2)(5,„.,,„+i] 
for TO = 1, 2, • • • ,p, whereas for the upper band. 



(5) 



1 yjm . 

— ;= cos a ^y sm a 



1 sin a 
sfn y 



pq 



E 



9V"- 
1 



^'A 



m'=p+2 y/im' -p){m' -p-l) 

X \9{m' - to) - [m - p - l)5m' ,ra] (6) 

for TO = p + 1, • • • , 71. We thus end up with the effective 



Hamiltonian in the diagonalized form 



dx 



2^ l^A 



[dx^Aix)f + VAKA[dxeAix)f 



(7) 

where v^^ = vi for = 1 , . . . , p — 2 and v^^ — u„ for 
v = p — 1^ . . . ,n — 2. The explicit formulae for the pa- 
rameters introduced in the transformation as well as the 
TL parameters are listed in the Appendix. Note that He 
{Ha) describes the Gaussian theory for charge excita- 
tions (inter-band excitations), for which the TL parame- 
ter Kc [Ka) characterizes the U(l) critical line. On the 
other hand, spin excitations in each band are described 
by -ff(T„, which still possess SU(p) (SU((7)) symmetry for 
the upper (lower) band. This fixes the value Ka^ = 1. 

2. 2 Expression for 1 /Ti at low temperatures 

We now derive the formula for the NMR relaxation 
rate, 1/Ti. Let us first write down the time-dependent 
spin-spin correlation function at finite temperatures, < 

S+{x,t)S-{Q,Q) >, 



< S+{x,t)S-{Q,Q) > 



En 



r ' 1 r / 

^ I sinh §^{x + v^t) sinh §^{x - v^t) 



2A- 



(8) 



with (3 = nkBT, where fi{= c. A, cr^) classifies massless 
modes in eq. (7). The phase factor e"'^ comes from the 
large momentum transfer across the Fermi points, where 
the momentum k{= 2kp'^\k'"p + A;^"') specifies such ex- 
citations. Here are conformal dimensions which are 
related to the scaling dimension At + A-. 

The dynamical spin susceptibility is expressed in terms 
of the above correlation functions, 

X±{k,u;)^-i J dtdxeit){[s+{x,t),s-{o,ow<^'^'-''-\ 

(9) 

The transverse spin correlation function, < S'^S'^ >, can 
be expressed in a similar manner. By computing eqs. (8) 
and (||) with the use of eqs. (5), (6) and (7), we obtain 
a generic formula for 1 /Ti at low temperatures. 



1 



= lim 



Ti c^-o n-'uj J 2tt 

^^i3,r,(A)r''»-i 



where 



r.(A) = 



n 



(10) 



(11) 
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Here A± is hyperfine form factor and B^, is a constant 
independent of the temperature. The corresponding crit- 
ical exponents ry^ = 2 ^ obtained as, 



T]M=2 h - <^ COS a -./-?; sin a ^ K, 




'7fc(o+fcW - 2 - 7 - - 



This completes the derivation of the fomiula for 1/Ti, 
which generalizes the results of Sachdevll3'' to multicom- 
ponent cases. It is to be noted here that we have explic- 
itly obtained the enhancement factor in (10), which 
has been neglected in ordinary literatures. We shall 
see that this quantity plays an important role for the 
NMR relaxation rate in ID systems, particularly when 
the system is located around the quantum phase transi- 
tion point around which the velocities are dramatically 
renormalized. Some explicit examples are shown in the 
following section. 

We have so far treated the model in the continuum 
limit within a weak coupling approach, so that it is not 
straightforward to evaluate 1/Ti as a function of the mi- 
croscopic model parameters. In the following section we 
evaluate the renormalized parameters and K^, by ex- 
ploiting the solvable model to deduce characteristic be- 
haviors in 1/Ti. 

§3. Exact Critical Properties of 1/Ti 

In this section, we exactly evaluate the renormalized 
parameters and to discuss critical properties of 
1 /Ti for the SU(7i) spin-orbital chain with band splitting. 
Also, the effect of hole-doping is discussed. 

3.1 Spin-orbital systems 

Let us start with the spin-orbital system in a insulating 
phase. In the limit U oo, the spin-orbital sector in the 
degenerate Hubbard model (1) is reduced to the SU(ri) 
antiferromagnetic Heisenberg model (each site is occu- 
pied by one electron), which is written down in terms of 



the fcrmion operators. 



(T.(T n .n 



(15) 

where only the single occupation of electrons is allowed at 
each site. The band splitting A due to, e.g. a crystalline 
field, forms two energetically separated bands, which are 
composed of the p lower and q upper bands including 
spin degrees of freedom, n = 2N = p + q. 

The exact solution of the model (15) was obtained 
by SutherlaHd, ^' and its properties have been clarified 
thus far.BB&lOEP Following standard procedures, we can 
evaluate the velocities of elementary excitations and the 
criticaJ—fixponents by applying finite-size scaling tech- 
niquesc^LlP to the excitation spectrumff 

In the insulating phase, charge excitations are frozen, 
while spin excitations are still gapless which possess 
SU(p) (SU(g)) symmetry for tkalower (upper) band. By 
applying the finite-size scalingEZP to the excitation spec- 
trum, we obtain the TL parameter as, 



-i 



pq 



A(p,9)' 



a 



0, 



(16) 



which generalizes the results of the SllMIjp&se.B' In the 
above expression, the dressed chargcE§E3'Ej'Ell' ^A(p,g) = 
CA(p,g) (-^p)j which features the critical line of the TL liq- 
uid, is given by the integral equation. 



CA(p,g)(Ap) = l+ / G(Ap - Ap)^A(p,9)(Ap)dAp,(17) 



where G(A) is defined as 

G(A) = Fp,p(A)+F,,,(A) 



1 

2^ 



dfc, 



^M,-(A) 



1 

2^ 



d/c. 



(18) 
(19) 



with X — exp(— |fc|/2). Here the cut-off parameters ztA^ 
are determined by-minimizing the free energy for a given 
band splitting A.Er 

Having determined the TL parameter Kj\, it is now 
straightforward to derive the critical exponents for var- 
ious correlation functions by exploiting the scaling rela- 
tions obtained in eqs.(12)-(14). Note that the TL param- 
eter Kc for the charge degree of freedom vanishes in the 
insulating system. Let us discuss the longitudinal spin 
susceptibility x_l, as an example. For excitations with 
momentum transfer 2kp'^^ within the lower and upper 



bands, we obtain the exponents 77 (, 
(13), 



from (12) and 



1 

1 - - 

P. 

1 

1 - - 

9. 



np 
nq 



(20) 



(21) 



On the other hand, for the momentum transfer kp^ ±kp^ 
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we use the relation (14), which yields 

1 1" 




(22) 



The scaling relations (20), (21) and (22), which are the 
generalizations of the ordinary TL relations to multicom- 
ponent cases, characterize critical properties of the de- 
generate model with band splitting. 

We note here that for the relaxation via excitations 
around the zero momentum transfer, the critical expo- 
nent is fixed to the canonical value, resulting in the spe- 
cific formula for this relaxation channel. 



1 



(23) 



which holds for any value of the band splitting, although 
the velocity v may be changed according to the splitting. 
Note that the quantity 1/u^, which originates from the 
enhancement factor in (10), is proportional to the 
square of the density of states for elementary excitations. 
Therefore, this expressicpL is regarded as a ID analogue 
of the Korringa relation.E3'' 



-1 1 1 r 



— ^=2kpC\'2kp^"'' ■ ' 




Fig. 1. Critical exponents of the NMR relaxation rate in the case 
of SU(7i) spin-orbital chain in which n-fold multiplet split into 
two set of n/2-fold states with the band splitting A: (a)n = 2 
(b)4, (c)6. Note that for the n = 2, 6 cases, A may be regarded 
as the magnetic field. The critical value Ac at which the upper 
band becomes massive is given by (a) 1.93, (b) 0.69 and (c) 0.36, 
respectively. 



Let us now focus on the specific case for which SU(n) 
multiplet splits into two SU(n/2) bands (n = 2iV). Note 
that the n = 4 case has been intensively studied thus 
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Fig. 2. The enhancement factor Fk in 1/Ti for the SU(4) chain 
as a function of the band splitting A. The critical field at which 
the upper band becomes massive is Ac = log 2. 



far.m^H'Bdi) In Fig. 1 we show the critical exponents as 
a function of the band splitting A. In the absence of the 
band splitting, A = 0, the TL parameter K^is given by 
— Ij which results in r]2kf. = 2(1 — l/n)n Hence the 
power-law temperature dependence becomes less singu- 
lar as the orbital degeneracy is increased. Although the 
critical exponents vary continuously when A is increased 
up to its critical value Ac, their values do not change so 
dramatically to induce the qualitative difference in the 
temperature dependence except for n — 2 (single orbital 
case). On the other hand, beyond Ac, the upper band 
becomes massive, being separated from the lower band. 
Thus for A > Ac, the low-energy physics can be de- 
scribed by the SU(7i/2) spin model. Correspondingly, 
the critical exponent shows a discontinuity at Ac. For 
example, r]^^(i) abruptly changes from 2(1 — l/p-fl/p^) to 

2(l-l/p) for the 811(0.) model at A = Ac. As claimed by 
Yamashita and Ueda,EP this discontinuity does not imply 
that the relaxation rate itself possesses such singularity, 
but there instead occurs the crossover behavior in 1/Ti. 
This point is discussed in more detail in the next section. 

We now observe how the enhancement factor F^ in 
(10) behaves as a function of the band splitting. The 
computed results for F^ are shown in Fig. 2. As A is in- 
creased, Fk is increased monotonically and is divergently 
enhanced near A = Ac. This characteristic behavior is 
caused by the renormalization of the velocities, which 
gives rise to the edge singularity in density of states. This 
may be regarded as a ID analog of the "mass enhance- 
ment" well known for the Fermi liquid in three dimen- 
sion. Note that Fq as well as F (o within the lower band 
still has the finite values even for A > Ac, whereas the 
others should vanish after divergently increased at Ac. 
As mentioned before, such a singularity is apparent, and 
should be replaced by the crossover behavior. We note 
that the above enhancement factor has been usually ne- 
glected to analyze the temperature dependence of 1/Ti, 
However, when the field- dependence of 1/Ti is studied 
at low temperatures, this should be properly taken into 
account, in particular around the critical point A = Ac. 
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3.2 Hole-doped systems 

We now turn to a hole-doped case. For a metallic 
system with finite holes, we consider the multicomponent 
t-J model, 



N 



H = -tJ2J2J2 (4.n' ^^+1-' + '^•^•) + Hj. (24) 

2=1 n'=l 

For the supersAopmetric case, t = J, the exact solution 
was obtained J3'eP which will be used here. Though this 
model is rather special in its appearance, some essential 
properties for a doped system can be described well, as 
is shown below. 

In Fig. 3, we show the critical exponents for the SU(n) 
symmetric case without the band-splitting A ^ (n = 
2,4,6). In this case, the critical exponent ri{— rj. 



2 2 ~ 
ri = 2-- + -K^ 
n n 



Kc = -i 

n 



2 



(25) 



where the dressed charge £,h(n) = C/i(n)(-^ri) is given by 



?h(n)(A„) = l+ / F„,„(A„ - A^)^5(„,„)(A^)dA',. (26) 

J -AO 

We note that the relation (25) generalizes the TL scal- 
ing relation to SU(ri) electron systems. The computed 
results are shown in Fig. 3. When the hole concentration 
is increased {fi is decreased), the critical exponent for 
the NMR relaxation rate, rj, is monotonically increased 
from 2(1 — 1/n + 1/n^) up to 2 (band bottom), because 
Kc — 1/n —f 1 in this case. We are now discussing a 
rather special solvable model, so that the global feature 
of the critical exponent should depend on the model. 
However, characteristic properties close to the insulating 
phase (/i ~ ^c) is expected to be rather general because 
in the small hole-doping region the charge sector is scaled 
to the strong coupling fixed point (hard-core bosons) ir- 
respective of the bare value of J. 
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Fig. 3. Critical exponents of the NMR relaxation rate for the 
SU(n) supersymmetric 4- Jmodel (n = 2, 4, 6) as a function of the 
chemical potential ^. At fj, = fic, the system becomes the Mott 
insulator; /ici ~ 1.39, fj.c2 ~ 0.44 and /ic3 ~ 0.22, respectively 
for n = 2, 4, 6. 



To see such properties for the critical exponents, we 
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Fig. 4. Critical exponents as a function of the band splitting A in 
a metallic system close to the Mott insulating phase: (a) SU(2), 
(b) SU(4) and (c) SU(6). 



now focus on the hole-doped case close to insulator (/i ^ 
/ic). In this limiting case, the TL parameters in (12), 
(13) and (14) are reduced to, 



Ka = 



y sma 




A(p,g) 



-i 



1 



1, 



0, (27) 



where in addition to the dressed charge (16), we have 
introduced another dressed charge S,B{p,q) = ^B{p,q){^n) 
which incorporates the interference between the charge 
and spin degrees of freedom. 



(A„) 



A2 



2,?(A„ - Ap)CB(p,9)(Ap)dAp.(28) 



In Fig. 4, the critical exponents in a metallic system 
with the small hole concentration are shown as a func- 
tion of the band splitting A {n = 2,4,6). As soon as 
holes are doped into the insulator, all the exponents are 
increased (compare Fig.l with Fig. 4), because the charge 
sector couples with the spin sector, and thus affects the 
relaxation process via "gauge interaction" due to Fermi 
statistics. For the limiting case, A Ac, the critical 
exponents approach 

„ 4 

V2kj.w = 2- 



n 

.2-1 
n 



16 
1 

2^ 



(29) 
(30) 
(31) 



Beyond this critical band-splitting, the NMR relaxation 
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rate is given by that for the SU(n/2) i-J model without 
band-splitting. Thus the critical exponent is discontinu- 
ously changed to smaller values at A = Ac for any case 
of n. 

We wish to mention that the coefficient is enhanced 
for small hole-doping, since the charge velocity Vc in 
eq.(ll) is very small in this region. However, it may 
not be easy to observe this enhancement due the renor- 
malized charge velocity, because the hole-doping itself 
increases the value of the exponents, thus having a ten- 
dency to hide singular behaviors at low temperatures. 

§4. Dynamical Spin Susceptibility 

We have shown that the critical exponents as well as 
the enhancement factor exhibit singularities around the 
critical point Ac. As already mentioned, such singular- 
ities may not be observed in NMR measurements, but 
should be replaced by the crossover behavior In order 
to clearly see how the crossover behavior emerges, we dis- 
cuss the dynamical spin susceptibility Imx(fc,Lj) in the 
low frequency regime. 

By exploiting the velocities and the TL parameters 
determined by the exact solution, we have computed 
lm.x{k,uj) at finite temperatures by using the expression 
(^). To be specific, we restrict ourselves to the SU(4) 
spin-orbital model with the band splitting. 

In Fig. 5, the spectral function Imx(fc, ct;)/^ is shown 
for A = at finite temperatures. For all the momenta 
shown in Fig. 5 (a), the broad maximum structure is 
found in the spectral function, indicating the presence of 
overdamped excitations. This is because we are dealing 
with the case of higher temperatures in (a), Tivq < ksT 
{v = 7r/2). On the other hand, if the temperature is de- 
creased with q being fixed, the sharp peak structure is 
developed around uj — Tivq in the spectral function (see 
(b) for which fivq 3> ksT) which indicates the presence 
of propagating spinous with the above dispersion rela- 
tion. These behaviors for the propagating mode axe anal- 
ogous to those observed for the SU(2) spin chainO Note 
that although there are three kinds of massless modes, 
all the velocities are the same for the SU(4) symmet- 
ric case with A = 0, so that we can observe only one 
spinon mode in the spectral function lmx{k , uj) / lu . The 
integration over q in the case of (b) determines the low- 
temperature power-law behavior of 1 /Ti discussed in the 
previous section. 

We next discuss the case with the finite band splitting 
close to its critical value Ac, where we have encountered 
the singularities in the NMR relaxation rate. In Fig. 6, 
we display the spectral function l'mx{k,uj)/u! by choos- 
ing two typical values of the temperature. There are 
three dominant low-energy excitations with different ve- 
locities around the momenta kp^ + k''p\= tt/2) 
and 2k p\ At higher temperatures shown in panels (a), 
(b) and (c), we can see the broadened maximum struc- 
ture at low frequencies, which is caused by overdamped 
spin excitations, as mentioned above. When the tem- 
perature is lowered, each excitation starts to show dif- 
ferent behaviors in the spectral function. We begin with 




0.02 0.04 0.06 0.08 0.1 



Fig. 5. Spectral function lmx{k,Lu)/LU of the SU(4) model (A = 
0) as a function of the frequency for several values of the mo- 
menta k{= it/2 + q): (a) higher temperature case (T = 0.005) 
and (b) lower temperature case (T = 0.0005). Note that the mo- 



(0 



2k 



mcntum q measures the difference from 7r/2{= 2k 
k^P + k^p^) around which elementary excitations become mass- 
less. 



the spectral function for the 2fc^'' mode in (d), and the 



k'"^' mode in (e). By comparing them with those 



,(0 , M 

r<,p -r l^p 

for the A = case, we find that the peak structure in the 
spectral weight is shifted to the extremely low frequency 
regime, and its peak-position is monotonically increased 
with the increase of the momentum q, reflecting the fact 
that the velocities of these excitations become very small 
around Ac. 

Recall here that the NMR relaxation rate (10) is de- 
termined by the integration of the spectral function over 
the momentum in the limit a> ^ 0. In the previous sec- 
tion, we have shown that 1/Ti is extremely enhanced 
by the velocity renormalization around Ac. However, 
we should keep in mind the following things. We are 
now considering the continuum limit of the original lat- 
tice model, so that when the system is close to the crit- 
ical point Ac, the energy range where such a treatment 
may be valid becomes pretty small. For instance, in the 
above case, the energy range where the dispersion rela- 
tion UJ 



hvuq holds should be in 



q ^ Itp' (the 

cutoff is given by k'p' ~ 0.15 for A = 0.66). Therefore 
the enhancement of 1/Ti can be seen in the very low fre- 
quency regime within the small cutoff parameter. Hence, 
even if such enhancement is naively expected to be ob- 
served at low temperatures, it may be obscured when the 
temperature is increased beyond the above cutoff energy. 

In contrast to the above two modes, the spectral func- 
tion for the 2fc^^ mode in Fig. 6(f) shows quite different 
behaviors; it develops two peaks of propagating spinous 
around 2kp\ whose dispersions are given by w ~ hvAQ 
and LU ~ hviq, respectively. When the system is close to 



NMR Relaxation Rate for ID Multicomponent Spin-Orbital Systems 



7 



A=0.66, T=0.005 



A=0.66, T=0.0005 



[xlO"] 4 

^ 3 

3 

s 1 



0.5 



[xlO- 



,^1 



3 6 

i 4 

& 2 

0, 



(a)2kF<"'+q 



q=0.005 
q=0.025 
q=0.045 



T — ' — I — ' — I — ' — r 



(b)kF'"'+kF"'+q 





0.02 0.04 0.06 0.08 0.1 0.02 0.04 0.06 0.08 0.1 



Fig. 6. Spectral function Imx{k,Lj)/Lj for the SU(4) model with 
the band splitting A = 0.66 (close to the critical value Ac ~ 
0.69). There are three massless modes, whose velocities are re- 
spectively given by ~ 2.87, Ca — 0.48 and Vu — 0.10. Note 
that the values of some solid lines have been scaled down by the 
factor of 1/5 or 1/2. 



the critical point, the interband excitations with the dis- 
persion v/^q becomes less important. On the other hand, 
in such cases, the excitation of w ~ hviq within the lower 
band becomes more dominant to control the relaxation 
process. In fact, beyond Ac, this excitation completely 
determines the l/Ti at low temperatures. 

From the above analyses of the dynamical spin suscep- 
tibility, we can see how the enhancement of 1/Ti in the 
vicinity of the critical field naturally crossovers to the re- 
gion where 1 /Ti is dominated by excitations in the lower 
band, when the energy or the temperature is raised. 

§5. Summary 

We have studied the NMR relaxation rate l/Ti for ID 
multicomponent spin-orbital systems. By generalizing 
the bosonization approach of Sachdev, we have obtained 
a generalized formula of 1 /Ti for multicomponent quan- 
tum systems with the band splitting. We have then ex- 
actly estimated the relaxation rate by using the SU(n) 
integrable model for the insulating as well as metallic 
phases. It has been found that the power-law temper- 
ature dependence of l/Ti becomes less singular as the 
orbital degeneracy is increased. This is also the case for 
the finite band splitting, as far as A is smaller than its 
critical value Ac. We have also pointed out that the re- 
laxation rate may be enhanced around Ac due to the 
dramatic renormalization of the velocities. This type of 
the enhancement may be common to ID quantum sys- 
tems when the spin/orbital gap is formed or collapsed by 



external fields. In particular, this effect becomes impor- 
tant to analyze the field- dependence of 1/Ti at low tem- 
peratures. It may be quite interesting to experimentally 
observe the enhancement in l/Ti, or the corresponding 
crossover behavior at very low temperatures in the NMR 
experiments. 
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Appendix 

The parameters for the canonical transformation in (5) 
and (6) are obtained as. 



tan2Q; — 



7 



I , . „ 1 

n 

qvi 



y' = 



+7 



- {pvi + qv^ 

'PVu 



6v 



pvi + qvu 

where vi and are the bare velocities, and the TL pa- 
rameters Ac and Aa are given by 

-1/2 

^ Z, 1 / /, — I H/,(l I > \ 



Ac 



Aa = 1 - 



TTVo 

2{n - l)aoC/ 

TTVo 



-1/2 



The above set of the parameters in the transformation 
gives the diagonalized Hamiltonian in the text, for which 
the renormalized velocities, Vc[a) ^-nd the TL parame- 



ters, A, 



c(A) 



read. 



Apq 
A2 



6v' 



1 

2 



7^ A^2 ^ ^ y(^) 



[Ac-2(1 + 27A^2) 



where T{a) (A^^ _ K^^)y/TTtai[2a. 
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